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In this paper we provide numerical evidence of the richer behavior of the connectivity 
degrees in heterogeneous preferential attachment networks in comparison to their homo- 
geneous counterparts. We analyze the degree distribution in the threshold model, a pref- 
erential attachment model where the affinity between node states biases the attachment 
probabilities of links. We show that the degree densities exhibit a power-law multiscaling 
which points to a signature of heterogeneity in preferential attachment networks. This 
translates into a power-law scaling in the degree distribution, whose exponent depends 
on the specific form of heterogeneity in the attachment mechanism. 
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1. Introduction 

A complex network is a set of nodes and links with a non-trivial topology (without ev- 
ident principles of design) d. Networks pervade all domains of science: examples can be 
found in the natural 121 social El, technological S and cultural El realms. In the current 
eff'ort to achieve a single coherent framework for complex systems, network theory has fo- 
cused on the underlying principles that govern their topology. The coupling architecture 
of a system is an important ingredient in its modeling because its topology aff'ects its 
collective dynamics and thus its function. Through a combination of empirical observa- 
tions, analytical insight and numerical sir nulat ion. complex networks have been subject 
to intensive scrutiny during the last years EEI. 

Dynamical network modelsElare discrete-time (generally stochastic) dynamical sys- 
tems that evolve networks by the iterated addition/subtraction of nodes/links. These 
regard network topology as an emergent property of its evolution, focusing on the mech- 
anisms that concur on such process. Among these mechanisms, preferential attachment 
enjoys a leading position in network literature. It was first discovered by Simon El in the 
1950s and applied to citation networks by PricelMIIl in the 1970s. Recently, the so-called 
preferential attachment model by Barabasi and Albert C21 has adopted this mechanism 
for the generation of undirected networks. The preferential attachment model prov ides 
a minimal account of sufficient mechanisms for the emergence of scale-free networks H^J 
Such networks are characterized by a connectivity distribution according to a power law, 
P(fc) = k which leads to a non-negligible presence of hubs (highly connected nodes). 

In this paper we present numerical results concerning the connectivity degrees in the 
threshold model. The threshold model belongs to a general class of heterogeneous prefer- 
ential attachment models where the affinity between node states biases the connectivity 
degrees in the attachment rule. First we review the basic facts about the preferential 
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attachment model. Next we introduce our framework for heterogeneous preferential at- 
tachment models, and we define the threshold model. Afterwards we briefly discuss the 
methodology followed in the numerical simulation and we discuss the results concerning 
the degree distribution and degree densities of the model, comparing them with the ana- 
lytical predictions made in the thermodynamic limit. We end this paper with concluding 
remarks about the significance and possible applications of the results. 



2. Preferential attachment 

The preferential attachment model by Barabasi and Albert li2J is a dynamical network 
model that prescribes two mechanisms in the evolution of a network: growth and pref- 
erential attachment. The process starts with a seed of arbitrary size, and a new node 
is added to the network at each step. Each newly added node has a fixed number m of 
links attached, which are connected to the already existing nodes following the so-called 
attachment rule. The rule states that the linking probability of a network node Vi is 
proportional to its connectivity degree ki, Il(vi) = ki/'^^ kj. The addition of nodes is 
iterated until a network with a desired size N is achieved. 

The properties of the preferential a ttachnient model have been widely studied, using 
both analytical and numerical methods 

I14I15I16I 

The networks generated exhibit power 
law deg ree distributions with 7 ~ 3, which becomes 7 = 3 in the thermodynamic limit 
N — > 00 li^ moderate clustering levels and short geodesic distances. The preferential at- 
tachment model is strictly topological as the node degrees are the only metric tha t drives 
network evolution. The model has been subject to a series of variations HZHSHS] that are 
still homogeneous in nature since the attachment rule doesn't take into account intrinsic 
properties of the nodes. Although network theory has led to a significant improvement 
in our understanding of complex systems, it has been argued that an augmentation of 
its framework is required in order to improve our modeling of complexity OEE 

More recently, heterogeneous networks have been proposed as a first step in address- 
ing the complication introduced by the infiuence of individual elements on the network 
structure. Weighted networks have been proposed i21)22 | 23| account for differences 
in the strengt h of the connecting links. Other proposed models include social netw ork 
models | 24 | 25 | make use of the so-called "social distance" ; a co mpetition model 

which incorporates a fitness measure; a loc al k nowledge model 12ZI which spatially lim- 
its the attachment of links; a metric model which make the attachment probabi lity 
proportional to a power of the Euclidean distance between nodes; a gas-like model 1221 
which interpr ets n odes as molecules that can exchange links upon collision; and topolog- 
ical automata where both node states and links destinations depend on the states of 
neighboring nodes in a previous iteration . 



3. Heterogeneous preferential attachment 

In order to enable a systematic analysis of the influence of heterogeneity in preferen- 
tial attachment networks we have introduced an extended formalism for heterogeneous 
models. The preferential attachment model can easily be generalized to heterogeneous 
networks by imposing a metric structure on the node states while preserving the orig- 
inal mechanisms of growth and preferential attachment. This provides an intermediate 
abstraction level between homogeneous networks and topological automata. 



3.1. Heterogeneous models with global affinity 

In this section we formally deflne a general class of heterogeneous models where node 
states bias the attachment probabilities of links in a generalized attachment rule. 



Definition 1. 

A heterogeneous preferential attachment model with global affinity Mi is a 3-tuple 
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(-R, p, a), where: 

(1) Each network node Vi is characterized by a state or attribute Xi on an arbitrary 
metric space R. 

(2) Each node state xi is randomly assigned following a distribution p over the state 
space R. 

(3) Each pair of nodes Vi and Vj has a nonnegative affinity a depending on the two state 



The class of heterogeneous preferential attachment models with global affinity HPAg is 
the set of all 3-tuples that satisfy the conditions in Definition 1: HPAg = {Mi : Mi = 
(R, p,cr)}. The formalism introduced defines the evolution of a network according to the 
following rules: 

(i) The node states describe intrinsic properties which are deemed constant in the 
timescale of evolution of the network. The network links are not characterized by 
any state or affinity. 

(ii) The growth process starts with a seed composed by Nq nodes and Lq links. The seed 
nodes Vi are assigned arbitrary states Xi £ R. 

(iii) A new node Va is added to the network at each iteration of the process. The newly 
added node Va has m links attached to it. The number m is common for all the added 
nodes and remains constant during the evolution of the network. The newly added node 
is randomly assigned a state Xa following the distribution p{x) over R. 

(iv) The m potential links attached to the newly added node Va are randomly connected 
to the network nodes following a probability distribution {n(s;i)} given by a generalized 
attachment rule, 

= (1) 

The visibility 7r of a node Vi in the attachment rule is given by the product of its 
connectivity degree ki and its affinity a with the newly added node Va, 

TT{Vi) = ki ■ a{Xi,Xa). (2) 

The function a measures the affinity of the interaction between nodes Vi and Va as a 
function of their states Xi and Xa- It thus can be seen in Eq. |2] that the affinity a for 
each potential interaction biases the connectivity degree ki of the candidate node in 
the attachment rule according to the states of the two nodes. The higher the affinity of 
a candidate node with the added node, the higher its probability Tl{vi) of acquiring a 
potential link. Steps (iii) and (iv) are iterated until a desired number of nodes has been 
added to the network. Obviously, only nodes with a strictly positive visibility can be 
regarded as receptive to the links attached to a newly added node. 



3.2. Threshold models 

The choice of the tuple {R, p, a) in the previous section determines the form of hetero- 
geneity in the attachment mechanism and thus the network topology. In order to apply 
numerical simulation we need to specify particular models within this general class. 
Next we define a family of heterogeneous preferential attachment models we refer as the 
threshold models. The threshold models are based on the assumption that the affinity 
between network nodes is proportional to the similarity of their intrinsic properties. In 
other words, the affinity a is inversely related to the distance between the node states 
as defined by the metric d on the state space R. In the spirit of the original preferential 
attachment model, the threshold models provide a minimal account of jointly sufficient 
mechanisms for the emergence of scale-free heterogeneous networks. 

As a first approximation we will consider the case where the transition of a is dis- 
continuous, this can be accounted by introducing a single parameter p that defines the 
threshold distance for the transition. Then the affinity is a reverse step function of the 
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Fig. 1. Affinity a in the threshold model (Eq. |4]l as a function of the distance d between states 
for a threshold /i = 0.4 and different width e values. The width values are: (a) e = 0, (b) e = 0.3. 



distance between the node states and can be regarded as a form of local attachment I2ZI 
adopting the expression: 

cr{Xi,Xa) = 1 - H{d{Xi,Xa) - At) = 1 - H^{d{Xi,Xa)), (3) 

where H{x) is the Heaviside function. The interaction threshold fi is shared by all the 
network nodes and governs the shape of the affinity in the preferential attachment rule, 
thus the form of heterogeneity in the network evolution. The shape of the affinity function 
a in this case is depicted in Fig. [T] (a) for a particular value of fi. 

Therefore, when the distance between the added node and the network node d{xi, Xa) 
is lower than the threshold fi, the receiving node is regarded as a valid candidate for 
the attachment and its affinity is maximum, cr = 1. In such a case the visibility of the 
network node is given by its connectivity degree, 7r(Di) = fc^. It should be noted that when 
the threshold /i is equal to the maximum possible distance, d„iax, there is obviously a 
maximum affinity between all the network nodes and the added node, irrespective of their 
states. In such circumstances the effective heterogeneity in the attachment mechanism 
is null, the visibility of all the network nodes is equal to their connectivity degrees and 
the dynamics recover the behavior of the homogeneous preferential attachment model 
proposed by Barabasi and Albert. 

On the other hand, when the distance d{xi,Xa) is larger than the threshold /i the 
affinity is inexistent, cr = 0. The lack of affinity between the nodes means that the 
visibility of the network node is zero from the standpoint of the newly added node, in 
this case their interaction is not possible since the attachment probability is Il(vi) = 
TT{vi) = 0. Last, when the threshold /i is small enough so that a newly added node 
cannot find a candidate node in the whole network, then it is considered that this node 
is rejected by the network and it is not accounted by the network size N. 

A more general case can be stated if we consider that the transition in the affinity 
(T takes place smoothly along a transition regime. For simplicity, let us approximate the 
inverse relationship between affinity and distance by means of a piece-wise linear map. 
This can be accounted by introducing an additional parameter e that defines the width 
of the transition regime. The shape of the affinity cr in this second case is depicted in 
Fig.[T](b). Thus we arrive at the following definition of a threshold model: 

Definition 2. 

A threshold model TMi S HPAg is a 3-tuple {R, p, cr) where cr : [0, dmax] [0, 1] is a 
nonnegative function defined for < ^ < dmax, < e < d 

max 

r 1 for < d < c^o 

a{xi,Xa) = (j{d{xi,Xa)) = < 1/2 + {fi - d)/€ for do < d < di (4) 
(0 for di < d < dmax, 
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where d is a metric on _R, dmax = sup^^ ^.^ g , a;2), dp = max(0, ^ — e/2) and 

di = mm(dmax,fi + e/2). 

It should bo noted that < cr < 1 , hence the affinity between nodes takes the highest 
value when the distance d < dg and the lowest value when d> di. The transition between 
these extreme affinity values takes place along a region defined by its middle point fi and 
its width e. Again, these two parameters are shared by all the network nodes and govern 
the shape of the affinity a in the preferential attachment rule. When the width e = 
the affinity becomes a step function as considered in the previous case. When e > the 
affinity cr can take values along the interval < cr < 1 and the transition in the visibility 
of the nodes between ki and is smooth. Likewise, rejected nodes are not accounted by 
the network size N. 

Henceforth we will focus on the unidimensional version of the threshold model over 
the real line, thus we will choose the unit interval R = [0, 1] as state space and the 
Euclidean distance d = \xi — Xa\ as metric of the space state. The affinity parameters fi 
and e of the model will therefore range in the interval [0, 1] . 

4. Connectivity degrees in the threshold model 

Next we present the results concerning the behavior of the connectivity degrees in the 
threshold model over the real line as a function of the affinity parameters fi and e. The 
connectivity degrees in the networks have been characterized by means of the degree 
distribution and the degree densities. These metrics have been analyzed by numerical 
simulation of ensembles of networks, each of them characterized by a particular pair of 
affinity parameters. The other model parameters are shared by all the network ensembles. 
The distributions have been computed by means of histograms of the relevant metrics 
in the different ensembles of networks. 

To implement the generalized attachment rule a vector of the visibilities of all in- 
cumbent nodes is calculated at each iteration according to Eq.[2] A list of the nodes with 
positive visibility is then compiled, and a candidate node in this list is randomly chosen 
for each new link with a probability given by the normalized visibility according to Eq. 
m Once a network node is selected and receives a new link at a given iteration, such node 
is removed from the list of candidate nodes during the remainder of the iteration. This 
follows from the fact that the network links are deemed as undirected and unweighted 
in the model, thus only one link can exist at most between a given pair of nodes. 

The simulations have been carried out for different thresholds fi ranging in the inter- 
val (0, 1], and for each threshold several widths have been considered e = {0, 0.1, 0.5}. 
The ensemble size is 10* networks, a sample size chosen to render an accurate estimation 
of the distributions around the cutoff points induced by finite-size effects. Each network 
has a final size of = 10* nodes, built from a random connected seed with A^o = 10 
nodes and Lq = 9 links, through the addition of nodes with m = 3 links and states 
uniformly distributed over the unit interval. In order to evaluate the robustness of the 
behavior of the model, the analysis has been repeated for different network sizes N and 
numbers of links m, without rendering qualitative discrepancies in the results. 

4.1. Degree distribution 

The degree distribution P(fc) of a network measures the probability for finding a node 
with degree k in the network. Alternatively, the product A^ ■ P{k) measures the average 
number of nodes in the network with a given connectivity degree k. As it has been 
stated before, it should be emphasized that the main asset of the preferential attachment 
model lies in its ability to generate networks with a scale-free degree distribution, that 
is, asymptotically following a power law P{k) ~ k~'^ . Next we study how robust is such 
scaling against changes in the form of the affinity cr of the threshold model, as defined 
by the parameters fi and e. 
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Fig. 2. Degree distribution A'^ ■ P{k) in the threshold model for different values of the threshold 
fj, and the width e. Rows (a), (b) and (c) correspond to /i = [0.5,1], [0.05,0.5] and (0,0.05], 
respectively. Columns (1) and (2) correspond to e = and 0.5, respectively. Within each graphic, 
a grayscale differentiates the values of the threshold /i. The other parameters are p{x) ^ U{0, 1), 
S = 10*, N = 10"', No = 10, Lo = 9 and m = 3. 
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Fig.[2]shows the results of the numerical simulation for the degree distribution P{k). 
Each curve depicted measures the average number of nodes A'^ • P{k) with degree k in 
a network generated with a particular vector of parameters. Alternatively, dividing the 
magnitude in the ordinate axis by the network size we obtain the average fraction of 
nodes with degree fc in a network. In the discussion of the results, let us focus first on 
the influence of /i when e = 0. This is the case when the affinity is a step function of the 
distance between the node states, cr = 1 — Hii{\xi — Xa\), and is represented in Fig. [2] 
Col. (1). Afterwards we will consider the effect of an increase in the transition width e. 

The curve /i = 1, depicted in Fig. |2] (l.a), illustrates the homogeneous case {/i = 
1) equivalent to the original preferential attachment model. This curve asymptotically 
follows a power law P{k) ^ k ' with an exponent 7 ~ 2.9. This value is equal to the 
exponent obtained by Barabasi and Albert in numerical simulations of the preferential 
attachment model (see for instance Fig. 21 in Ref.El), The slight discrepancy in the value 
of 7 with regards to the analytical prediction HSl is a finite-size effect that decreases the 
effective value of the exponent as the network size is reduced. Likewise, for larger degrees 
the distribution undergoes a crossover to an exponential decay, which is another finite- 
size effect due to the limited number of nodes in the simulated networks. The step-like 
plateaux are also a finite-size effect due to the relatively poorer statistics of hubs (that 
is, nodes with very high connectivity). 

When the threshold value fi is modified, the numerical results in Fig. |2] (l.a)-(l.c) 
show that the power-law scaling in the degree distribution P{k) is robust against changes 
in the heterogeneity of the attachment mechanism. The results also evidence that the 
variation in the scaling exponent 7 induced by changes in the form of heterogeneity 
is slight. These two traits exhibit a fair agreement of the numerical results with the 
behavior of the a naly tical solution for the stationary degree distribution of the model. In 
a previous paper we showed that in the thermodynamic limit the degree distribution 
P{k) of a heterogeneous model M = (R, p, cr) in the HPAg class adopts for k > m the 
expression: 

=f I U ^^^] 

Jr \^J-J-+, + 2 J wm + 2 
where w{x) is a normalized fitness factor, 

wix) = = , - - . ■ , with«)(x) = / cr{x,y)p{y)dy. 

w Jj^w(x)p(x)dx Jr 

For a convenient comparison of the analytical and numerical results, Fig. |3] depicts 
the stationary degree distribution obtained from Eq. [5]for the threshold model, next to 
the degree distribution obtained by numerical simulation, when e = 0. A first evident 
trait in the numerical results shown in Fig.[2](l.a) and (l.b) is the existence of a change 
of regime in the influence of fi observed around fi = 0.5. When fi > 0.5, Fig. [2] (l.a) 
evidences that a decrease in the threshold p. yields a lower scaling exponent (7 ~ 2.9 for 
p = I, while 7 ~ 2.6 for p = 0.5) and a higher cutoff degree in the distribution P{k), 
which points towards a higher presence of hubs. 

The behavior exhibited by the exponent 7 of the numerical distribution P{k) as a 
function of a decreasing p agrees with the behavior of the stationary distribution, which 
was analytically explained as a consequence of the introduction of density components 
f{k,x) with a slower decay, which dominate the asymptotical behavior of P{k). Alter- 
natively, this behavior can be explained through the increasing inequality of the affinity 
regions of the network nodes over the state space R, which translates into an increase 
in the visibility of nodes with moderate states and a decrease in the visibility of nodes 
with peripheric states. As the threshold p approaches 0.5, the differences in the visibility 
of the nodes grow, so that the network nodes with largest affinity regions face a lower 
competition in the acquisition of links and attain higher degrees. 

It should be noted that as the threshold p approaches 0.5, the exponent 7 in the 
analytical solution exhibits a slightly better agreement with the numerical results. This 
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can be attributed to an interplay between finite-size effects and the inhomogeneity of the 
density components of P(k). On the one hand, the decrease of the network size tends to 
magnify finite-size effects that yield a reduction in the exponent 7. On the other hand, 
the decrease of the network size tends to reduce the relative advange accrued by the 
most visible nodes, and thus yield an increase in the same exponent 7. Therefore, the 
inequality in the affinity regions brings a better agreement between the stationary and 
numerical distributions. 

When the threshold ^ < 0.5, Fig. [2](l-b) evidences that a decrease in still yields a 
lower exponent 7 (7 ~ 2.5 for = 0.05) as it was observed when ^ > 0.5. Furthermore, 
a decrease in ^ also yields a reduction in the cutoff degree to the exponential regime, 
which points to a lower presence of hubs in the networks. The lower presence of hubs in 
the simulated networks is in agreement with the behavior of the stationary distribution, 
however the behavior of the exponent 7 shows a discrepancy with regards to the analytical 
predictions. The regime reversal exhibited by the presence of hubs in the networks was 
analytically explained as a consequence of the progressive removal of density components 
f{k, x) with slower decay. 

Alternatively, it should be noted that a reduction of fj, below 0.5 yields a decrease 
in the inequality of the affinity regions of the network nodes over the state space R, 
which yields a weakening of the advantage accrued by the most connected nodes and a 
subsequent reversal in the cutoff degree. Concurrently, a decrease of fi also yields a grow- 
ing localization of the attachment mechanism over the state space R and a progressive 
degradation of the interaction domains of the newly added nodes. The fragmentation in 
the network dynamics translates into a magnification of the finite-size effects, a weak- 
ening of the positive feedback mechanism in the degree dynamics, and thus a further 
left shift in the cutoff degree. The growing finite-size effects also yield a reduction in the 
exponent 7, which dominates over the increase in 7 produced by the increasing equality 
in the balance of competition between the network nodes. 

Finally, when the threshold fi <^ 1, Fig. |2](l-c) shows that the numerical distribu- 
tion P{k) continues the trend exhibited for larger thresholds in Fig. [2] (l-b). The scaling 
exponent 7 decreases slightly and the cutoff degree shifts to the left. In the limit when 
fi approaches zero, P{k) converges to a power-law with an exponent 7 ~ 2.2. It should 
be emphasized that in the last case the finite-size effects introduce noticeable differ- 
ences in the behavior of the numerical distribution against the stationary distribution, 
as evidenced by the lower scaling exponent and the reduced cutoff degree. In the ther- 
modynamic limit N — > 00, a reduction of the threshold from fj, = 0.5 to — ► yields an 
increase in the scaling exponent from 7 ~ 2.6 to 7 3, as shown in Fig. |3] (l-b). 

The variable influence of the finite-size effects on the degree distribution P{k) as a 
function of the threshold value fj, can be more clearly appreciated in Fig. 13] which illus- 
trates the results of the numerical simulation for different network sizes and two threshold 
values. Clearly, larger network sizes yield higher cutoff degrees for the exponential regime 
irrespective of the threshold fj, adopted. On the other hand, the influence of the network 
size A'^ on the scaling exponent of P{k) largely depends on the threshold value fi, which 
governs the degree of fragmentation of the preferential attachment on the state space. 
When fi takes moderate or high values, as in the case of Fig.|4](l), an increase in the size 
A'^ yields a slight increase in the exponent 7. As the threshold fj, decreases and approaches 
zero, which is the case of Fig. |4](2), flnite-size effects become more important and an 
increase in A'^ yields a progressively larger increase in the exponent 7. 

To sum up, the growing inequality in the affinity regions of the network nodes is 
the main factor driving the observed phenomena when fi > 0.5, while the dilution of 
the former inequality and the fragmentation of the attachment mechanism are the main 
factors driving the observed phenomena when fi < 0.5. These factors alter the balance 
of competition for links in the preferential attachment mechanism and the subsequent 
presence of hubs in the networks. It should be emphasized that only the last factor has 
a finite-size origin, thus the agreement between the thermodynamic calculations and the 
numerical simulations (with reasonable network sizes) deteriorates as the heterogeneity 
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k k 

Fig. 3. Comparison of the analytical solution and the numerical simulation of the degree dis- 
tribution P{k) in the threshold model. Rows (a) and (b) correspond to ^ = [0.5,1] and (0,0.5], 
respectively. Column (1) corresponds to the analytical solution (Eq. [Sjl with p{x) ~ U{0, 1), e = 
and m = 3. Column (2) corresponds to the numerical estimation with p{x) ~ U{0, 1), e = 0, 
S = 10", N = 10", A^o = 10, Lo = 9 and m = 3. 

level increases. These results prompts us to consider the heterogeneity in the attachment 
mechanism as an influential ingredient in the hub composition of real networks. 

Next we scrutinize the influence of the transition width e in the numerical degree 
distribution P{k). When e > the affinity is no longer a step function of the distance 
d between node states, as illustrated in Fig.[l](b). Despite the qualitative change in the 
form of (T, the simulation results furnish evidence that when e is small enough (e.g. e = 
0.01) the behavior of the distribution P{k) is barely distinguishable from that obtained 
when e = 0, irrespective of the value of the threshold /i (results not shown). Therefore, 
the influence of the threshold fi on the degree distribution P{k) is robust against the 
functional form of cr, in the sense that slight changes in e do not introduce noticeable 
differences in the behavior of the distribution. 

As the width e takes higher values, the numerical results prove that the power- 
law scaling in the degree distribution P{k) is preserved, as predicted by the analytical 
study. The results also show that as e increases the sensitivity of P{k) to changes in 
the threshold /i progressively decreases. When the threshold fj, = 1 the distribution is 
virtually the same for all widths e, however as /i is reduced from 1 to 0.5 the variation 
experienced by P{k) is smaller as e increases, as illustrated in Fig. [2] (2. a). As /i is 
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Fig. 4. Degree distribution N ■ P{k) in the threshold model for different network sizes A'^ and 
two values of the threshold fi. Columns (1) and (2) correspond to = 0.5 and 10"'^, respectively. 
Within each graphic, a grayscale differentiates the network sizes N . The other parameters are 
e = 0, p(x) ~ (7(0, 1), S = 10*, Nq = 10, Lq = 9 und m = 3. 



reduced below 0.5 the variation experienced by P{k) is even smaller for the lower half of 
the threshold interval, as illustrated in Fig. [2] (2. b). This is particularly evident for very 
small values of the threshold fi, as represented in Fig. |2] (2.c). 

The influence of the width e on the stationary distribution P{k) of the threshold 
model was analytically explained by the introduction of quadratic terms in the normal- 
ized fitness w, which tend to dampen the dependence of the lowest scaling exponents 7 
in the density components of P{k) as a function of the threshold fi for higher e values. 
Alternatively, this phenomenon may be explained by the shape of the affinity function 
a pictured in Fig. [T] As the width e increases, the transition regime in n becomes slower 
and this performs a sort of averaging on the interaction domains of the added nodes: 
more nodes are included in the domains, however the visibility of some is downgraded by 
the transition regime. This averaging effect dampens the dependence of the distribution 
P{k) as a function of the threshold p., in particular for small thresholds fj, when the shape 
of the cr function becomes dominated by the transition regime. 



4.2. Degree density 

The degree density /(fc, x) of a heterogeneous network measures the probability density 
for finding a node with degree k and state x in the network. Assuming an uncountable 
set of states over the space R, as in the case of the threshold model studied in this paper, 
then f^^ f{k, x)dx measures the probability for finding a node with degree k and state 
a;i < X < X2. Our analytical studies of heterogeneous preferential attachment models 
l2iJ suggest that the presence of heterogeneity in the attachment mechanism leaves a 
signature in the topology of the networks that can be empirically checked. This signature 
is a power-law multiscaling in the degree densities according to exponents spanning a 
continuous interval around the homogeneous exponent 7 = 3. The stationary density 
components f{k, x) adopt the expression for k > ni: 

f(k X) = ^^/^ B{k,l + 2M 

m + 2/w B{Tn,l + 2/w)' ^' 



where B is the Legendre Beta function, 



B{y,z)=f t»-^{l-ty-^dtlovy,z>0. 
Jo 



(8) 
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Fig. 5. Comparison of the analytical solution and the numerical simulation of the degree densities 
f(k, x) in the threshold model for different node states x. Column (1) corresponds to the analytical 
solution (Eq. ^ with p{x) ~ U{0, 1), fi = 0.5, e = and m = 3. Column (2) corresponds to the 
numerical approximation by means of the partial degree distributions P{k\x) with p{x) ~ U{0, 1), 
fj, = 0.5, € = 0, m = 3 and = 0.01. Within each graphic a grayscale differentiates the states x 
that characterize each subset of network nodes. The other parameters are as in Fig. |3] 

This means that the asymptotic behavior of f{k,x) when fe — > oo is as 

/(fc.x) ~B(fc, l + 2/«;) (9) 

so that the scaling exponents adopt the form: 

7(x) = 1 + 2/w)(x). (10) 

Next we present the results of the simulation of the degree densities that validate 
this prediction in the case of the threshold model and provide a numerical evidence of 
the existence of the multiscaling phenomenon. In order to compute the degree densities 
f{k,x) we have estimated partial degree distributions P(k\x) over subsets of network 
nodes with states in an interval on certain xq values of width , x g (xq — Tx-.xq + r^). 
In the limit when the width Tx approaches zero, the partial distributions P{k\x) converge 
to the degree densities f{k,x). The procedure has been carried out for different xq 
values ranging in the [0.5, 1] interval with a width = 0.01, small enough to obtain 
an accurate approximation to the degree densities. Due to symmetry considerations the 
[0, 0.5) interval has not been scrutinized. 
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Fig. [5] shows the stationary degree densities /(fc, x) obtained from Eq. [7] for the 
threshold model, next to the partial degree distributions P{k\x) obtained by numerical 
simulation, when e = 0. Each curve depicted in Fig. |5] col. (2) measures the probability 
for finding a node with with connectivity degree k and state x satisfying the inequality 
\x — xo\ < Tx in a. network of the ensemble. It can be noted that the partial distributions 
exhibit a power-law scaling with different exponents that depend on the state that char- 
acterizes the subset of nodes. The scaling exponent is highest (7 ~ 3.9) for the subset 
of nodes with states around the most peripheric sample state, xq = 0.99. The exponent 
decreases as the sample state becomes more moderate, and is lowest (7 ~ 2.4) for the 
subset of nodes with states around the central state x = 0.5. The numerical results 
exhibit a good agreement with the analytical predictions in Fig. [5] col. (1) and provide 
evidence of the multiscaling phenomenon observed in the stationary degree densities. 

5. Conclusions 

In summary, in this paper we have studied by numerical simulation the scaling phe- 
nomena exhibited by the connectivity degrees in the threshold model. We have shown 
that the power-law scaling in the degree distribution is a robust trait in the presence of 
heterogeneity, however significant details such as the scaling exponent or the presence 
of hubs in the network depend on the form of such heterogeneity in the attachment 
mechanism. We have also shown that the results of the numerical simulations exhibit a 
fair agreement with the stationary distributions of heterogeneous preferential attachment 
networks analytically obtained in the thermodynamic limit. The existence of finite-size 
effects introduces slight discrepancies in the numerical results that only increase as the 
threshold value becomes very small. 

The robustness evidenced by the scaling phenomenon in the degree distribution 
points out that heterogeneity provides a more general mechanism for the emergence of 
power-law distributions, with regards to the original preferential attachment mechanism 
suggested by Barabasi and Albert. The values of the scaling exponents (2.2 — 2.9 in the 
numerical simulations, 2.5 — 3.0 in the analjrtical calculations) also fit nicely within the 
range of values empirically observed in real networks. We suggest that this fact may help 
explain the ubiquity of real scale-free networks that otherwise differ in other topological 
aspects. The presence of heterogeneity thus may constitute a source of differentiation 
behind the behavior exhibited by less universal metrics. 

Finally, we have shown that the presence of heterogeneity induces a multiscaling 
phenomenon in the degree densities of preferential attachment networks. Such multiscal- 
ing is not found in their homogeneous counterparts and hence represents a signature 
of heterogeneity in the topology of preferential attachment networks. We have shown 
that this phenomenon can be experimentally observed by computing the partial degree 
distributions on subsets of network nodes with states in small neighborhoods of the state 
space. This result provides an experimental procedure to check the presence of hetero- 
geneity in the evolution of real networks when this process is suspected of having arisen 
from some form of preferential attachment. 
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